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Abstract. Is a spontaneous transition from a state x to a state y allowed by thermodynamics? 
Such a question arises often in chemical thermodynamics and kinetics. We ask the more formal 
question: is there a continuous path between these states, along which the conservation laws hold, the 
concentrations remain non-negative and the relevant thermodynamic potential G (Gibbs energy, for 
example) monotonically decreases? The obvious necessary condition, G(x) > G(y), is not sufficient, 
and we construct the necessary and sufficient conditions. For example, it is impossible to overstep 
the equilibrium in 1-dimensional (ID) systems (with n components and n— 1 conservation laws). The 
system cannot come from a state x to a state y if they are on the opposite sides of the equilibrium even 
if G(x) > G(y). We find the general multidimensional analogue of this ID rule and constructively 
solve the problem of the thermodynamically admissible transitions. 

We study dynamical systems, which are given in a positively invariant convex polyhedron D and 
have a convex Lyapunov function G. An admissible path is a continuous curve in D along which 
G does not increase. For x,y £ D, x ^ y (x precedes y) if there exists an admissible path from x 
to y and x ~ y if x £ y and y £ x. The tree of G in D is a quotient space D/ ~. We provide an 
algorithm for the construction of this tree. In this algorithm, the restriction of G onto the 1-skeleton 
of D (the union of edges) is used. The problem of existence of admissible paths between states is 
solved constructively. The regions attainable by the admissible paths are described. 

Key words. Lyapunov function, convex polyhedron, attainability, tree of function, entropy, free 
energy 
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1. Introduction. 

1.1. Motivation, ideas and a simple example. "Applied dynamical sys- 
tems" are models of real systems. The available information about the real systems 
is incomplete and uncertainties of various types are encountered in the modeling. Of- 
ten, we view them as errors: errors in the model structure, errors in coefficients, in 
the state observation and many others. Nevertheless, there is an order in this world 
of errors: some information is more reliable, we trust in some structures more and 
even respect them as laws. Some other data are less reliable. There is an hierarchy 
of reliability, our knowledge and beliefs (described, for example by R. Peierls [53] for 
model making in physics). Extracting as many consequences from the more reliable 
data either without or before use of the less reliable information is a task which arises 
naturally. 

In our paper, we study dynamical systems with a strictly convex Lyapunov func- 
tion G defined in a positively invariant convex polyhedron D. For them, we analyze 
the admissible paths, along which G decreases monotonically, and find the states that 
are attainable from the given initial state along the admissible paths. The main area 
of applications of these systems is chemical kinetics and thermodynamics. The mo- 
tivation of our research comes from the hierarchy of reliability of the information in 
these applications. 

Let us discuss the motivation in more detail. In chemical kinetics, we can rank 
the information in the following way. First of all, the list of reagents and conservation 
laws should be known. Let the reagents be A±, A2, . . . , A n . The non- negative real 
variable Ni > 0, the amount of A{ in the mixture, is defined for each reagent, and TV 
is the vector of composition with coordinates Ni. The conservation laws are presented 
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by the linear balance equations: 

n 

bi(N) = ^^Nj = const (i = l,...,m). (1.1) 
i=i 

We assume that the linear functions b{(N) (i = 1, . . . , m) are linearly independent. 

The list of the components together with the balance conditions (|1.1|) is the first 
part of the information about the kinetic model. This determines the space of states, 
the polyhedron D defined by the balance equations ([1.1)) and the positivity inequalities 
Ni > 0. This is the background of kinetic models and any further development is 
less reliable. The polyhedron D is assumed to be bounded. This means that there 
exist such coefficients that the linear combination Xibi(N) has strictly positive 
coefficients: \{a\ > for all j = 1, . . . , n. 

The thermodynamic functions provide us with the second level of information 
about the kinetics. Thermodynamic potentials, such as the entropy, energy and free 
energy are known much better than the reaction rates and, at the same time, they 
give us some information about the dynamics. For example, the entropy increases in 
isolated systems. The Gibbs free energy decreases in closed isothermal systems under 
constant pressure, and the Helmholtz free energy decreases under constant volume 
and temperature. Of course, knowledge of the Lyapunov functions gives us some 
inequalities for vector fields of the systems' velocity but the values of these vector 
fields remain unknown. If there are some external fluxes of energy or non-equilibrium 
substances then the thermodynamic potentials are not Lyapunov functions and the 
systems do not relax to the thermodynamic equilibrium. Nevertheless, the inequality 
of positivity of the entropy production persists and this gives us useful information 
even about the open systems. Some examples are given in [26) [28], 

The next, third part of the information about kinetics is the reaction mechanism. 
It is presented in the form of the stoichiometric equations of the elementary reactions: 

^apiAi^Y^PpiAi, (1.2) 

i i 

where p = 1, . . . , m is the reaction number and the stoichiometric coefficients a P i, j3 P i 
(i = 1, . . . , n) are nonnegative integers. 

A stoichiometric vector j p of the reaction ([1.2)) is a n-dimensional vector with 
coordinates 7^ = f3 p i — a P i, that is, 'gain minus loss' in the pth elementary reaction. 

The concentration of A{ is an intensive variable q = Ni/V, where V > is the 
volume. The vector c = N/V with coordinates q is the vector of concentrations. A 
non- negative intensive quantity, r p , the reaction rate, corresponds to each reaction 
(|1.2p . The kinetic equations in the absence of external fluxes are 

p 

If the volume is not constant then equations for concentrations include V and have 
different form. 

For perfect systems and not so fast reactions the reaction rates are functions of 
concentrations and temperature given by the mass action law and by the generalized 
Arrhenius equation. A special relation between the kinetic constants is given by the 
principle of detailed balance: For each value of temperature T there exists a positive 
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equilibrium point where each reaction (|1.2p is equilibrated with its reverse reaction. 
This principle was introduced for collisions by Boltzmann in 1872 [10]. Wegscheider 
introduced this principle for chemical kinetics in 1901 [67 . Einstein in 1916 used it in 
the background for his quantum theory of emission and absorption of radiation [T7] . 
Later, it was used by Onsager in his famous work [51 . For a recent review see [50] . 

At the third level of reliability of information, we select the list of components 
and the balance conditions, find the thermodynamic potential, guess the reaction 
mechanism, accept the principle of detailed balance and believe that we know the 
kinetic law of elementary reactions. However, we still do not know the reaction rate 
constants. 

Finally, at the fourth level of available information, we find the reaction rate 
constants and can analyze and solve the kinetic equations (|1.3|) or their extended 
version with the inclusion of external fluxes. 

Of course, this ranking of the available information is conventional, to a cer- 
tain degree. For example, some reaction rate constants may be known even better 
than the list of intermediate reagents. Nevertheless, this hierarchy of the information 
availability, list of components - thermodynamic functions - reaction mechanism - 
reaction rate constants, reflects the real process of modelling and the stairs of available 
information about a reaction kinetic system. 

It seems very attractive to study the consequences of the information of each 
level separately. These consequences can be also organized 'stairwise'. We have the 
hierarchy of questions: how to find the consequences for the dynamics (i) from the list 
of components, (ii) from this list of components plus the thermodynamic functions of 
the mixture, and (iii) from the additional information about the reaction mechanism. 

The answer to the first question is the description of the balance polyhedron D. 
The balance equations ([l.ip together with the positivity conditions Ni > should 
be supplemented by the description of all the faces. For each face, some Ni = and 
we have to specify which Ni have zero value. The list of the corresponding indices 
z, for which Ni = on the face, / = {ii, . . . , i^}, fully characterizes the face. This 
problem of double description of the convex polyhedra [49] [14] [21] is well known in 
linear programming. 

The list of vertices [6] and edges with the corresponding indices is necessary 
for the thermodynamic analysis. This is the 1-skeleton of D. Algorithms for the 
construction of the 1-skeletons of balance polyhedra as functions of the balance values 
were described in detail in 1980 [26]. The related problem of double description for 
convex cones is very important for the pathway analysis in systems biology [58[ [22] . 

In this work, we use the 1-skeleton of D, but the main focus is on the second step, 
i.e. on the consequences of the given thermodynamic potentials. For closed systems 
under classical conditions, these potentials are the Lyapunov functions for the kinetic 
equations. For example, for perfect systems we assume the mass action law. If the 
equilibrium concentrations c* are given, the system is closed and both temperature 
and volume are constant then the function 



is the Lyapunov function; it should not increase in time. The function G is pro- 
portional to the free energy F = RTG + const (for detailed information about the 
Lyapunov functions for kinetic equations under classical conditions see the textbook 
[68] or the recent paper [33]). 




(1.4) 
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G(c*) A 3 A 2 



Fig. 1.1. The balance simplex (a), the levels of the Lyapunov function (b) and the thermody- 
namic tree (c) for the simple system of three components, Ai, A2, A3. Algorithm for finding a vertex 
v>zc (d). 

If we know the Lyapunov function G then we have the necessary conditions for 
the possibility of transition from the vector of concentrations c to c' during the non- 
stationary reaction: G(c) > G(c') because the inequality G(c(to)) > G{c{t$-\-t)) holds 
for any time t > 0. 

The inequality G(c) > G(c') is necessary if we are to reach d from the initial 
state c by a thermodynamically admissible path, but it is not sufficient because in 
addition to this inequality there are some other necessary conditions. The simplest 
and most famous of them is: if D is one-dimensional (a segment) then the equilibrium 
c* divides this segment into two parts and both c(to) and c(to + t) (t > 0) are always 
on the same side of the equilibrium. 

In ID systems the overstepping of the equilibrium is forbidden. It is impossible 
to overstep a point in dimension one, but it is possible to circumvent a point in 
higher dimensions. Nevertheless, in any dimension the inequality G(c) > G(c') is not 
sufficient if we are to reach c' from the initial state c along an admissible path. Some 
additional restrictions remain in the general case as well. A two-dimensional example 
is presented in Fig. 11.11 Let us consider the mixture of three components, ^-1,2,3 with 
the only conservation law c\ + C2 + C3 = b (we take for illustration 6 = 1) and the 
equidistribution in equilibrium c\ = c\ = c\ — 1/3. The balance polyhedron is the 
triangle (Fig. II. Ik ). In Fig. II. lb the level sets of 

3 

G = 5^Ci(ln(3ci)-l) 
i=i 

are presented. This function achieves its minimum at equilibrium, G(c*) = — 1. On 
the edges, the function G achieves its conditional minimum, go, in the middles, and 
go = ln(3/2) — 1. G reaches its maximal value, g max = In 3 — 1, at the vertices. 

If G(c*) < g < go then the level set G(c) = g is connected. If go < g < Vax then 
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the corresponding level set G(c) = g consists of three components (Fig. II. lb ). The 
critical value is g = go. The critical level G(c) = go consists of three arcs. Each arc 
connects two middles of the edges and divides D in two sets. One of them is convex 
and includes two vertices, the other includes the remaining vertex. 

A thermodynamically admissible path is a continuous curve along which G does 
not increase. Therefore, such a path cannot intersect these arcs 'from inside', i.e. 
from values G(c) < go to bigger values, G(c) > go- For example, if an admissible path 
starts from the state with 100% of A2, then it cannot intersect the arc that separates 
the vertex with 100% A\ from two other vertices. Therefore, any vertex cannot be 
reached from another one and if we start from 100% of A2 then the reaction cannot 
overcome the threshold ~77.3% of Ai, that is the maximum of c\ on the corresponding 
arc (Fig. II. lb ). This is an example of the 2D analogue of the ID prohibition of 
overstepping of equilibrium. 

For x,y G D, x >^y (x precedes y) if there exists a thermodynamically admissible 
path from x to and x ~ y if x £3 y and y £3 x. The equivalence classes with 
respect to x ~ y in D are the connected components of the level sets G(c) = g. 
The quotient space T = D/ ~ is the space of these connected components. For the 
canonical projection we use the standard notation ir : D — ^ T. This is the tree of the 
connected components of the level sets of G. (Here "tree" stands for a one dimensional 
continuum, a sort of dendrites [13], and not for a tree in the sense of the graph theory.) 

If x ~ y then G{x) = G(y). Therefore, we can define the function G on the tree: 
G(tt(c)) = G(c). It is convenient to draw this tree on the plane with the vertical 
coordinate g = G(x) (Fig. II. lb ). The equilibrium c* corresponds to a root of this 
tree, 7r(c*). If G(c*) < g < go then the level set G(c) = g corresponds to one point on 
the tree. The level G(c) = go corresponds to the branching point, and each connected 
component of the level sets G(c) = g with go < g < g m ax corresponds to a separate 
point on the tree. The terminal points ("leaves" with g > go) of the tree correspond 
to the vertices of D. 

An ordered segment [x, y] or [y, x] (x ^3 y) on the tree T consists of such points z 
that x £3 z )ZV- A continuous curve ip : [0, 1] —> D is an admissible path if and only if 
its image nocp : [0, 1] — ^ T is a path that goes monotonically down in the coordinate g. 
Such a monotonic path in T from a point x to the root is just a segment [x, 7r(c*)]. On 
this segment, each point y is unambiguously characterized by g = G{y). Therefore, 
if for c G D we know the value G(c) and a vertex v £3 c, then we can unambiguously 
describe the image of c on the tree: 7r(c) is the point on the segment [tt(v), 7r(c*)] with 
the given value of G, g = G(c). 

We can find a vertex v £3 c by a chain of central projections: the first step is the 
central projection of c onto the border of D with center c*. The result is the point c' 
on a face (in Fig. II. Id this is the point c' on an edge). The second step is the central 
projection of the point c' onto the border of the face with the center at the partial 
equilibrium c* r (that is, the minimizer of G on the face) and so on (Fig. II. lb ). If the 
projection on a face is the partial equilibrium then for any vertices v of the face v £3 c. 
In particular, if the face is a vertex v then v ^3 c. For a simple example presented in 
Fig. II. Id this is the vertex A\. 

In this paper, we extend these ideas and observations to any dynamical system, 
which is given in a positively invariant convex polyhedron and has there a strictly 
convex Lyapunov function. The class of chemical kinetic equations for closed systems 
provides us standard and practically important examples of the systems of this class. 
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1.2. A bit of history. It seems attractive to use an attainable region instead 
of the single trajectory in situations with incomplete information or with information 
with different levels of reliability Such situations are typical in many areas of science 
and engineering. For example, the theory for the continuous-time Markov chain is 
presented in [2, 27 and for the discrete-time Markov chains in [3]. 

Perhaps, the first celebrated example of this approach was developed in biological 
kinetics. In 1936, A.N. Kolmogorov [40] studied the dynamics of interacting popula- 
tions of prey (x) and predator (y) in the general form: 

x = xS(x,y), y = yW(x,y) 

under monotonicity conditions: dS(x,y) / 'dy < 0, dW(x,y)/dy < 0. The zero iso- 
clines, given by equations S(x, y) = or W(x, y) = 0, are graphs of two functions 
y(x). These isoclines divide the phase space into compartments with curvilinear 
borders. The geometry of the intersection of the zero isoclines, together with some 
monotonicity conditions, contain important information about the system dynamics 
that we can find [3D] without exact knowledge of the kinetic equations. This approach 
to population dynamics was applied to various problems [45j [7] . The impact of this 
work on population dynamics was analyzed in the review [62] . 

In 1964, Horn proposed to analyze the attainable regions for chemical reactors 
[36] . This approach became popular in chemical engineering. It was applied to the 
optimization of steady flow reactors [23] , to batch reactor optimization without knowl- 
edge of detailed kinetics [19], and for optimization of the reactor structure [34J. An 
analysis of attainable regions is recognized as a special geometric approach to reactor 
optimization [TS] and as a crucially important part of the new paradigm of chemical 
engineering [35] . 

Many particular applications were developed, from polymerization [63] to particle 
breakage in a ball mill [47] and hydraulic systems [28] . Mathematical methods for the 
study of attainable regions vary from Pontryagin's maximum principle [46J to linear 
programming [38], the Shrink- Wrap algorithm [43], and convex analysis. In 1979 it 
was demonstrated how to utilize the knowledge about partial equilibria of elementary 
processes to construct the attainable regions [24] . The attainable regions significantly 
depend on the reaction mechanism and it is possible to use them for the discrimination 
of mechanisms [29] . 

Thermodynamic data are more robust than the reaction mechanism. Hence, there 
are two types of attainable regions. The first is the thermodynamic one, which use 
the linear restrictions and the thermodynamic functions [25 . The second is generated 
by thermodynamics and stoichiometric equations of elementary steps (but without 
reaction rates) [2H[3T]. R. Shinnar and other authors [61] rediscovered this approach. 
There was even an open discussion about priority [9 . 

Some particular classes of kinetic systems have rich families of the Lyapunov 
functions. Krambeck [41] studied attainable regions for linear systems and the l\ 
Lyapunov norm instead of the entropy. Already simple examples demonstrate that 
the sets of distributions which are accessible from a given initial distribution by linear 
kinetic systems (Markov processes) with a given equilibrium are, in general, non- 
convex polytopes [241 [271 ED] • The geometric approach to attainability was developed 
for all the thermodynamic potentials and for open systems as well [26J. Partial results 
for chemical kinetics and some other engineering systems are summarized in [681 HE] • 

The tree of the level set components for differentiate functions was introduces in 
the middle of the 20 century by Adelson- Velskii and Kronrod [TJ H2] and Reeb [56] . 
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Sometimes these trees are called the Reeb trees [20] but from the historical point of 
view it may be better to call them the Adelson-Velskii - Kronrod - Reeb (or AKR) 
trees. These trees were essentially used by Kolmogorov and Arnold [4 in solution of 
the Hilbert's superposition problem (the ideas, their relations to dynamical systems 
and role in the Arnold's scientific life are discussed in his lecture [5]). 

The general Reeb graph can be defined for any topological space X and real 
function / on it. It is the quotient space of X by the equivalence relation defined 
by x ~ y holds if and only if f(x) = f(y) and x, y are in the same connected 
component of f~ 1 (f(x)). Of course, this "graph" is again not a discrete object from 
the graph theory but a topological space. It has application in differential topology 
(Morse theory [55]). in topological shape analysis and visualization [20l[39], in data 
analysis [64] and in asymptotic analysis of fluid dynamics [44, 59] . The books [20) [39] 
include many illustration of the Reeb graphs. The efficient mesh-based methods for 
the computation of the graphs of level set components are developed for general scalar 
fields on 2- and 3-dimensional manifolds [T6] . 

Some time ago the tree of entropy in the balance polyhedra was rediscovered as an 
adequate tool for representation of the attainable regions in chemical thermodynamics 
[25j[26]. It was applied to analysis of various real systems [37j[69]. Nevertheless, some 
of the mathematical backgrounds of this approach were delayed in development and 
publications. Now, the thermodynamically attainable regions are in extensive use in 
chemical engineering and beyond [SIig^OSEllE 

63, 69]. In this paper we aim to provide the complete mathematical background for the 
analysis of the thermodynamically attainable regions. For this purpose, we construct 
the trees of strictly convex functions in a convex polyhedron. This problem allows a 
general meshless solution in higher dimensions because topological and geometrical 
simplicity (the domain D is a convex polyhedron and the function G is strictly convex 
in D). In this paper, we present this solution in detail. 

1.3. The problem of attainability and its solution. Let us formulate pre- 
cisely the problem of attainability and its solution before the exposition of all technical 
details and proofs. Our results are applicable to any dynamical system that obeys a 
continuous strictly convex Lyapunov function in a positively invariant convex polyhe- 
dron. The situations with uncertainty, when the specific dynamical system is not given 
with an appropriate accuracy but the Lyapunov function is known, give a natural area 
of application of these results. 

Here and below, D is a convex polyhedron in R n , Do consists of the vertices of 
D, D\ is the union of the closed edges of D, that is, the 1-skeleton of D, and D\ is 
the graph whose vertices correspond to the vertices of D and edges correspond to the 
edges of D, (the graph of the 1-skeleton) of D. We use the same notations for vertices 
and edges of D and D\. 

Let a real continuous function G be given in D. We assume that G is strictly 
convex in D [57 . Let x* be the minimizer of G in D and let g* — G(x*) be the 
corresponding minimal value. 

The level set S g = {x G D\ G(x) = g} is closed and the sublevel set U g = {x G 
D | G(x) < g} is open in D (i.e. it is the intersection of an open subset of R n with 

Let us transform D\ into a labeled graph. Each vertex v G Do is labeled by the 
value j v = G(v) and each edge e = [v, w] C D\ is labeled by the minimal value of G on 
the segment [v,w] C D, g e = min^^] G(x). The vertices and edges of D\ are labeled 
by the same numbers as the correspondent vertices and edges of D\. By definition, 
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the graph Di \ U g consists of the vertices and edges of D\, whose labels 7 > g. 

The graph D\ \ U g depends on g but this is a piecewise constant dependence. It 
changes only at g — 7, where 7 are some of the labels of the graph D\. Therefore, 
it is not necessary to find this graph and to analyze connectivity in it for each value 

G(y)=g. 

Definition 1.1. A continuous path <p[0, 1] — » D is admissible if the function 
G((p(x)) does not increase on [0, 1]. For x,y G D, x £3 y (x precedes y) if there exists 
an admissible path (p[0, 1] —> D with (p(0) — x and tp(l) — y; x ~ y if x £3 y and 

The relation "£3" is transitive. It is a preorder on D. The relation is an 
equivalence. 

Definition 1.2. The tree of G in D is the quotient space T= D/ ~. 

The equivalence classes of ~ are the path-connected components of the level sets 
S g . For the natural projection of D on T we use the notation ir : D — » T. We denote 
by tt~ 1 (z) C D the set of preimages of z G T. The preorder on D transforms 
into a partial order on T: ir(x) £3 7r(?/) if and only if x ^ y. We call T also the 
thermodynamic tree keeping in mind the thermodynamic applications. The "tree" T 
is a ID continuum. We have to distinguish this continuum from trees in the graph- 
theoretic sense which have the same graphical representation but are discrete objects. 
In Sec. 13.21 ( "Coordinates on the thermodynamic tree") we describe the tree structure 
of this continuum. It includes the root, the edges, the branching points and leaves 
but the edges are represented as the real line segments. 

Definition 1.3. Let x,y G T, x ^3 y. An ordered segment [x,y] (or [y,x]) 
consists of such points z G T that x £3 z ^3 y. 

In Sec. [3] we prove that any ordered segment [x, y] (x ^ y) in T is homeomorphic 
to [0, 1]. A continuous curve (p : [0, 1] —> D is an admissible path if and only if its 
image ir o cp : [0, 1] — >• T is monotonic in the partial order on T. Such a monotonic 
path in T from x to y (x ^3 y) is just a path along a segment [x, y\. Each point z on 
this segment is unambiguously characterized by the value of G(z). 

We also use the notation [x, y] for the usual closed segments in R n with ends x,y: 
[x,y] = {Xx + (1 — \)y I A G [0, 1]}. The degenerated segment [x,x] is just a point 
{x}. The segments without one end are (x,y] and [x,y) and (x,y) is the segment in 
R n without both ends. 

The attainability problem: Let x,y G D and G(x) > G(y). Is y attainable from x 
by an admissible path? 

The solution of the attainability problem can be found in several steps: 

1. Find two vertices of D, v x and v y , that precede x and correspondingly. 
Such vertices always exist. There may be several such vertices in D. We can 
use any of them. 

2. Construct the graph Di \ Uc( y ) by deletion from D\ all the elements with the 
labels 7 < G(y). 

3. y is attainable from x by an admissible path if and only if v x and v y are 
connected in the graph D\ \ Ug( v )- 

So, to check the existence of an admissible path from x to y we should check the 
inequality G{x) > G(y) (the necessary condition) then go up in G values and find the 
vertices, v x and v yi that precede x and y, correspondingly (such vertices always exist). 
Then we should go down in G values to G(y) and check whether the vertices v x and v y 
are connected in the graph D\ \ Uc{ y )- The classical problem of determining whether 
two vertices in a graph are connected may be solved by many search algorithms 
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[52 , 50 , for example, by the elementary breadth-first or depth-first search algorithms. 

The procedure "find a vertex v x G Do that precedes x G D" can be implemented 
as follows: 

1. If x = x* then any vertex v G Do precedes x. 

2. If x ^ x* then consider the ray r x = {x* + A(x — x*) | A > 0}. The intersection 
r x DD is & closed segment [x*, x'\. We call x' the central projection of x onto 
the border of D with the center x' ^ x. 

3. The central projection x r always belongs to an interior of a face D r of D, 
< dimD' < dimD. If dim!) 7 > then set x := x' D := L>', x* := 
argmin{G(z) | z G and go to step 1. 

4. If dimD' = then it is a vertex v ^ x we are looking for. 

The dimension of the face decreases at each step, hence, after not more than dim D — 1 
steps we will definitely obtain the desired vertex. A simple example is presented in 
Fig. fTTH . _ 

The information about all connected components of D\ \ U g for all values of 
g is summarized in the tree of G in D, T (Definition II. 2p . The tree T can be 
described as follows (Theorem I3.3j) : it is the space of pairs (g,M), where g G 
[mino G(x), max/) G(x)] and M is a connected component of Di\U g , with the partial 
order relation: (g, M) £ (g* M f ) if g > g' and M C M' . For x,y e D, x ^ y if and 
only if tt(x) £ 7r(y). 

The tree T may be constructed gradually, by descending from the maximal value 
of G, g = g max (Sec. 13.3)) . At g = ^ max , the graph D\ \ U g consists of the isolated 
vertices with the labels 7 = g max (generically, this is one vertex). Going down in 
g, we add to D\ \ U g the elements, vertices and edges, in descending order of their 
labels. After adding each element we record the changes in the connected components 
of2?i\C^. 

For each point z G T, z = (g, M), its preimage in D, 7r _1 (^, M), may be described 
by the equation G(x) = g supplemented by a set of linear inequalities. Computation- 
ally, these linear inequalities can be produced by a convex hull operation from a finite 
set. This finite set is described explicitly in Sec. 13.41 

For each point z = (g,M) the set of all z' = (g',M') attainable by admissible 
paths from z has a simple description, g' < g, M' D M . 

The tree of G in D provides a workbench for the analysis of various questions 
about admissible paths. It allows us to reduce the n-dimensional problems in D to 
some auxiliaryjjuestions about such ID or even discrete objects as the tree T and the 
labeled graph D\. For example, we use the thermodynamic tree to solve the following 
problem of attainable sets: For a given x G D describe the set of all y ^ x by a system 
of inequalities. For this purpose, we find the image of x in T, 7r(x), then define the 
set of all points attainable by admissible paths from tt(x) in T and, finally, describe 
the preimage of this set in D by the system of inequalities (Sec. 13.4)) . 

1.4. The structure of the paper. In Sec.[2j we present several auxiliary propo- 
sitions from convex geometry. We constructively describe the result of the cutting of 
a convex polyhedron D by a convex set U : The description of the connected compo- 
nents of D \ U is reduced to the analysis of the ID continuum D\ \ U, where D\ is 
the 1-skeleton of D. 

In Sec.[3j we construct the tree of level set components of a strictly convex function 
G in the convex polyhedron D and study the properties of this tree. The main result of 
this section is the algorithm for construction of this tree (Sec. 13.3)) . This construction 
is applied to the description of the attainable sets in Sec. 13.41 These sections include 
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some practical recipes and it is possible to read them independently, immediately after 
Introduction. Several examples of the thermodynamic trees for chemical systems are 
presented in Sec. 2J 

2. Cutting of a polyhedron D by a convex set U. 

2.1. Connected components of D \ U and of Pi \ U. Let P be a convex 
polyhedron in W 1 . We use the notations: AS(D) is the minimal linear manifold that 
includes D; d = dim AS (D) = dimD is the dimension of P; ri(D) is the interior of 
D in Aff(P); rd(D) is the border of D in AS(D). 

For P, Q C W 1 the Minkowski sum is P + Q = {x + y \ x G P, y G Q}. The convex 
hull (conv) and the conic hull (cone) of a set V C W 1 are: 



convi 




cone(V) = < ^ XiVi 



q > 0, vi,...,v q G V, Ai,...A g > 0, 



For a set P C W 1 the following two statements are equivalent (the Minkowski-Weyl 
theorem) : 

1. For some real (finite) matrix A and real vector 6, D = {x G M n |Ax < b} ; 

2. There are finite sets of vectors {^i, . . . , v q } C W 1 and {n, . . . r p } C M n such 
that 

P> = conv{^i, . . . v g } + cone{r*i, . . . ,r p } . (2.1) 

Every polyhedron has two representations, of type (1) and (2), known as (halfspace) 
i^-representation and (vertex) V-representation, respectively. We systematically use 
both these representations. Most of the polyhedra in our paper are bounded, therefore, 
for them only the convex envelope of vertices is used in the V-representation (|2.ip . 
The fc-skeleton of D, D^, is the union of the closed /c-dimensional faces of D: 

DoCDiC...cD d = D. 

Do consists of the vertices jof D and D\ is a one-dimensional continuum embedded in 
R n . We use the notation D\ for the graph whose vertices correspond to the vertices 
of D and edges correspond to the edges of D, and call this graph the graph of the 
1-skeleton of D. 

Let U be a convex subset of R n (it may be a non-closed set). We use Uo for the 
set of vertices of D that belong to C/, Uo = U D P^ and U\ for the set of the edges of 
D that have non-empty intersection with U. By default, we consider the closed faces 
of P, hence, the intersection of an edge with U either includes some internal points 
of the edge or consists from one_of its ends. We use the same notation U\ for the set 
of the corresponding edges of D\. 

A set W C P C W 1 is a path- connected component of P if it is its maximal path- 
connected subset. In this section, we aim to describe the path-connected components 
of D \ U. In particular, we prove that these components include the same sets of 
vertices as the connected components of the graph Di \ U. This graph is produced 
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from Di by deletion of all the vertices that belong to Uq and all the edges that belong 
to Ut. 

Lemma 2.1. Let x G D \ U. Then there exists such a vertex v G Dq that the 
closed segment [v,x] does not intersect U: [v,x] C D \ U. 

Proof. Let us assume the contrary: for every vertex v G Z^o there exists such 
X v G (0, 1] that x + X v (v — x) G U. The convex polyhedron D is the convex hull 
of its vertices. Therefore, x = ^2 veDo k v v fc> r some numbers n v > 0, v G Do, 

Let 

A — ^ v 
It is easy to check that ^2 veD S v = 1 and 

veD 

According to (|2.2p , x belongs to the convex hull of the finite set {x + X v (v — x) \ v G 
Do } C U. U is convex, therefore, x G U but this contradicts to the condition x £ U. 
Therefore, our assumption is wrong and there exists at least one v G Dq such that 

[v,x] n u = 0. □ 

So, if a point from the convex polyhedron D does not belong to a convex set U 
then it may be connected to at least one vertex of D by a segment that does not 
intersect U. Let us demonstrate now that if two vertices of D may be connected in D 
by a continuous path that does not intersect U then these vertices can be connected 
in D\ by a path that is a sequence of edges D, which do not intersect U. 

Lemma 2.2. Let v,v' G D , v,v' £ U. Suppose that <p : [0,1] -> (D\U) is 
a continuous path, <p(0) = v and <p(l) = v' . Then there exists such a sequence of 
vertices {^o, . . . ,vi} C (D\U) that any two successive vertices, Vi,Vi+\, are connected 
by an edge e M+i C (D 1 \ U). 

Proof. Let us, first, prove the statement: the vertices v,v' belong to one path- 
connected component of D\U if and only if they belong to one path- connected com- 
ponent of D\\U . 

Let us iteratively transform the path (p. On the kth iteration we construct a path 
that connects v and v' in D^-k \ U, where d = dimD and k = 1, . . . , d — 1. We start 
from a transformation of path in a face of D. 

Let S C Dj be a closed j-dimensional face of D, j > 2 and let i(j : [0, 1] — » (Dj\U) 
be a continuous path, -0(0) = v, = v' and i/j([Q, 1]) H U = 0. We will transform 
ip into a continuous path ' [0, 1] —> (Dj \ U) with the following properties: (i) 
^(0) = v, ^(1) = v', (h) ^([0, 1]) H U = 0, (hi) ^([0, 1])\SC ^([0, 1]) \ S and 
(iv) -05 ([0,1]) fl ri(S) = 0. The properties (i) and (ii) are the same as for -0, the 
property (hi) means that all the points of ^s([0, 1]) outside S belong also to -0([0, 1]) 
(no new points appear outside S) and the property (iv) means that there are no points 
of ^sQO, 1]) in ri(S). To construct this ips we consider two cases: 

1. U H ri(S) ^ 0, i.e. there exists y° G U H ri(S); 

2. Unri(S) =0. 

In the first case, let us project any i/j(t) G ri(S) onto rd(S) from the center y°. Let 
y G 5, y y°. There exists such a X(y) > 1 that y° + X(y)(y — y°) G rd(S). This 
function A(y) is continuous in S \ {y }. The function X(y) can be expressed through 
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the Minkowski gauge functional [32] defined for a set K and a point x: 



p K (x) = inf{r > | x G rK}; X(y) = (p( D - yo )(y ~ Vo)) 



-l 



Let us define for any y G ri(S), y ^ y° a projection ns(y) = y° + ^(y)(y — 2/°)- This 
projection is continuous in 5 \ {y }, and 71-5(2/) = y if y E rd(S). It can be extended 
as a continuous function onto whole Dj\{y }: 



The center y° G U. Because of the convexity of U, if y ^ J7 then y° -\-X(y — y°) £ U 
for any A > 1. Therefore, the path ^s(£) = 7rs (?/>(£)) does not intersect J7 and satisfies 
all the requirements (i)-(iv). 

Let us consider the second case, U D ri(S') = 0. There are the moments of the 
first entrance of ip{t) in S and the last going of this path out of S: 



< r\ < T2 < 1. Let 2/ 1 = V ; ( r i) an d ?/ 2 = ^( r 2)- If 2/ 1 = 2/ 2 then we can just delete 
the loop between n and T2 from the path ip{r) and get the path that does not enter 
ri(S). So, let y 1 ^ y 2 . 

These points belong to rd(S). Let y s G ri(S) be an arbitrary point in the relative 
interior of S which does not belong to the segment [y 1 , y 2 ] (dim S > 2). The segments 
[y 1 ^ 8 ] an d [y 2 iy S ] do not intersect U because the following reasons: U D S C rd(S) 
(may be empty), neither nor belong to U, and all other points of the 3- vertex 
polygonal chain [y 1 ^ 5 \y 2 ] belong to ri(S). L 

Let P(y 1 1 y s \y 2 ) be a plane that includes the chain [y 1 , y s , y 2 ]. The intersection 
S fl P(y x , y s \y 2 ) is a convex polygon. The convex set U D 5 H P(y x , y s \y 2 ) belongs 
to the border of this polygon. Therefore, it belongs to one side of it (Fig. 12. ip (may 
be empty) because convexity of the polygon and of the set U. The couple of points 
y 1 ^ 2 cut the border of the polygon in two connected broken lines. At least one of 
them does not intersect U (Fig. 12. ip . Let us substitute ip on the interval [ti,T2] by 
this broken line. The new path does not intersect ri(S). Let us use for this new path 
the notation ips(t)- The path ips does not intersect ri(S) and J7, and all the points 
on them outside S are the points on the path ip for the same values of the argument 



So, for any closed face S C D with dim S = j > 2 and a continuous path 
ip : [0, 1] — » (Dj \ U) that connects the vertices v, v' of D (ip(0) = v, ip{l) = v') we 
construct a continuous path ips : [0, 1] — >• {Dj \ U) that connects the same vertices, 
does not intersect ri(S) and takes no new values outside 5, "0s ([0, C ^([0, 

Let us order the faces S C D with dim S > 2 in such a way that dim Si > dim Sj 
for i < j: D = So, 5i, . . . , Si. Let us start from a given path <p : [0,1] —> D\U that 
connects the vertices v and v' and let us apply sequentially the described procedure: 



By the construction, this path does not intersect any relative interior ri(Sk) (k = 
0, 1, . . . ,£). Therefore, the image of 6 belongs to D u : [0, 1] -» (Di \ U). It can be 
transformed into a simple path in D\ \ U by deletion of all loops (if they exist). This 
simple path (without self-intersections) is just the sequence of edges we are looking 




Ti = min{r | ip(r) GS}, T2 = max{r | ip{r) G S} , 



T . 



= {---(((<Ps o )sM...)s t . 



for. □ 
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Fig. 2.1. Intersection of a face S with the plane P(y 1 ,y s ,y 2 ) when U Dri(S) = (Lemma \2.2\ 
case 2). In this intersection, U n S C rd(S) belongs to one side of the polygon (the bold segment). 
The dashed lines with arrows represent the 3-vertex polygonal chain [y 1 , y s , y 2 ] . There exists a path 
from y 1 to y 2 along the boundary of the polygon. In Fig., this is the polygonal chain that follows 
the solid lines with arrows. 

Lemmas 12. 1[ 12.21 allow us to describe the connected components of the d-dimen- 
sional set D\U through the connected components of the one-dimensional continuum 
Dt\U. 

Proposition 2.3. Let Wi, . . . , W q be all the path- connected components of D\U . 
Then Wi D Z^o 7^ for all i = 1, . . . , q, the continuum Di\U has q path- connected 
components and Wi D D\ are these components. 

Proof. Due to Lemma 12. If each path-connected component of D \ U includes 
at least one vertex of D. According to Lemma 12. 2\ if two vertices of D belong to 
one path-connected component of D \ U then they belong to one path-connected 
component of D\ \ U . The reverse statement is obvious, because D\ C D and a 
continuous path in D\ is a continuous path in D. □ 

We can study connected components of a simpler, discrete object, the graph D\. 
The path-connected components of D \ U correspond to the connected components of 
the graph D\ \ U . (This graph is produced from D\ by deletion all the vertices that 
belong to /To and all the edges that belong to U\). 

Proposition 2.4. Let Wi, . . . , W q be all the path- connected components of D\U . 
Then the graph D\\U has exactly q connected components and each set Wi D Do is 
the set of the vertices of D of one connected component of Di\U. 

Proof. Indeed, every path between vertices in Di includes a path that connects 
these vertices and is the sequence of edges. (To prove this statement we just have to 
delete all loops in a given path.) Therefore, the vertices vi,V2 belong to one connected 
component of D\ \ U if and only if they belong to one path-connected component of 
D\ \ U . The rest of the proof follows from Proposition 12.31 □ 

We proved that the path-connected components of D \ U are in one-to-one corre- 
spondence with the components of the graph D\\U (the correspondent components 
have the same sets of vertices). In applications, we will meet the following problem. 
Let a point x G D \ U be given. Find the path-connected component of D \ U which 
includes this point. There are two basic ways to find this component. Assume that we 
know the connected components of D\ \ U. First, we can examine the segments [x, v] 
for all vertices v of D. At least one of them does not intersect U (Lemma 12. ip . Let it 
be [x,fo]. We can find the connected component Di\U that contains ^o- The point 
x belongs to the correspondent path-connected component of D \ U. This approach 
exploits the ^-description of the polyhedron D. The work necessary for this method 
is proportional to the number of vertices of D. 
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Another method is based on projection on the faces of D. Let x G ri(D). We can 
take any point y° G D\U and find the unique Ai > 1 such that x 1 = yo + \i(x — y°) G 
rd(D). Let x 1 G ri(Si), where Si is a face of D. If Si H U = then we can take 
any vertex vo G Si and find the connected component Di\U that contains vq. This 
component gives us the answer. If Si H U ^ then we can take any y 1 G Si H C/ and 
find the unique A2 > 1 such that x 2 = 2/ 1 -f A2(x 1 — y 1 ) G rd(S). This x 2 belongs to the 
relative boundary of the face Si. If x 2 is not a vertex then it belongs to the relative 
interior of some face S2, dimS2 > and we have to continue. At each iteration, 
the dimension of faces decreases. After d = dimD iterations at most we will get the 
vertex v we are looking for (see also Fig. 11.1)) and find the connected component of 
Di\U which gives us the answer. Here we exploit the i7-description of D. 

2.2. Description of the connected components of D \ U by inequalities. 

Let Wi, . . . , W q be the path-connected components of D \ U. 

Proposition 2.5. For any set of indices I c {1, . . . , q} the set 

is convex. 

Proof. Let y 1 ,y 2 G Kj. We have to prove that [y 1 ^ 2 ] C Kj. Five different 
situations are possible: 

1. y\y 2 eU; 

2. y 1 GU,y 2 GW z ,ieI; 

3. y\y 2 eWi,ieI, [y\y 2 ]nU = tt; 

4. y\y 2 eWi,ieI, [y\y 2 ]nU^9; 

5. y 1 G Wi, y 2 G Wj, ij G I, i^j. 

We will systematically use two simple facts: (i) the convexity of U implies that its 
intersection with any segment is a segment and (ii) if x 1 G Wi and x 2 G D \ Wi then 
the segment [x 1 ,^ 2 ] intersects U because Wi is a path-connected component of U. 
In case 1, [y 1 ^ 2 ] C U C K because convexity U. 

In case 2, there exists such a point y 3 G (y 1 ^ 2 ) that [y 1 ^ 3 ) CUD [y 1 ^ 2 ] C 
[y 1 ^ 3 ]. The segment (y 3 ,y 2 ] cannot include any point x G D \ Wi because it does 
not include any point from U. Therefore, in this case (y 3 ,y 2 ] C Wi C K and y 3 G K 
because it belongs either to U or to Wi. 

In case 3, [y 1 , y 2 ] C Wi C K because Wi is a path-connected component of D \ U 
and [y 1 ^ 2 ] n U = 0. 

In case 4, [2/ 1 ,?/ 2 ] H J7 is a segment L with the ends x x ,x 2 . It may be [x x ,x 2 ] 

< x 1 < x 2 < y 2 ), (x 1 ^ 2 } (y 1 < x 1 < x 2 < y 2 ), [x x ,x 2 ) (y 1 < x 1 < x 2 < y 2 ), 
or (x l ,x 2 ) (y 1 < x 1 < x 2 < y 2 ). This segment cuts [y 1 ^ 2 ] in three segments: 
[y 1 ^ 2 ] = L\ ULUL2, Li includes 2/ 1 and L2 includes y 2 . Therefore, L\ C Wi, L C J7 
and L2 C Wi because is a path-connected component of D \ U and U is convex. 
So,[y\y 2 ]cK. 

In case 5, [y 1 , y 2 ] HU is also a segment L with the ends x 1 , x 2 . It may be [x 1 , x 2 ] 
(y 1 < x 1 < x 2 < y 2 ), (x 1 ^ 2 } (y 1 < x 1 < x 2 < y 2 ), [x x ,x 2 ) (y 1 < x 1 < x 2 < y 2 ), 
or (x 1 ,^ 2 ) (y 1 < x 1 < x 2 < y 2 ). This segment cuts [y 1 ^ 2 ] in three segments: 
[y 1 ^ 2 ] = Li ULUL2, I/i includes y 1 and L2 includes ?/ 2 . Therefore, L\ <zWi, L <zU 
and L 2 C Wj because Wi j are path-connected components of D \ U and U is convex. 
So, [y\y 2 ] CK. □ 
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Typically, the set U is represented by a set of inequalities, for example, G(x) < g. 
It may be useful to represent the path-connected components of D\ U by inequalities. 
For this purpose, let us first construct a convex polyhedron Q C U with the same 
number of path-connected components in D \ Q : Vi, . . . , V q and with inclusons Wi C 
Vi. We will construct Q as a convex hull of a finite set. Let us select the edges e of D 
which intersect U but the intersection eC\U does not include vertices of D. For every 
such edge we select one point x e G efl U. The set of these points is Q\. By definition, 

Q = conv(/7 UQi) . (2.3) 

Q is convex, hence, we can apply all the previous results about the components of 
D \ U to the components of D \ Q. 

Lemma 2.6. The set Uo U Qi is the set of vertices of Q. 

Proof A point x G Uo^JQi is not a vertex of Q = conv(L r oUQi) if and only if it is a 
convex combination of other points from this set: there exist such xi, . . . , x^ G UoUQi 
and Ai , . . . , Afe > that Xi ^ x for all i = 1, . . . , k and 

k k 

Aj = i , XiXi = x . 

2=1 i=l 

If x G C/o then this is impossible because x is a vertex of D and C/q U Qi C If 
x e Qi then it belongs to the relative interior of an edge of D and, hence, may be a 
convex combination of points D from this edge only. By construction, Uo U Qi may 
include only one internal point from an edge and in this case does not include a vertex 
from this edge. Therefore, all the points from Qi are vertices of Q. □ 

Lemma 2.7. The set D\Q has q path- connected components Vi, . . . , V q that may 
be enumerated in such a way that Wi C Vi and Wi = Vi \ U. 

Proof To prove this statement about the path-connected components, let us 
mention that Q and j7 include the^same vertices of D, the set [To, and cut the same 
edges of D. Graphs Di\Q and Di\U coincide. Q C U because of the convexity of 
U and definition of Q. To finalize the proof, we can apply Proposition 12.41 □ 

Proposition 2.8. Let I be any set of indices from {1, . . . , q}. 

q u (jj = c ° nv u u (y ( d ° n f ) ) ( 2 - 4 ) 

Proof On the left hand side of (|2.4p we see the union of Q with the connected 
components Vi (i G I). On the right hand side there is a convex envelope of a finite 
set. This finite set consists of the vertices of Q, (Uo U Qi) and the vertices of D that 
belong to Vi (i G I). Let us denote by Ri the right hand side of (|2.4p and by Lj the 
left hand side of (j2^4|) . 

L/ is convex due to Proposition [23] applied to Q and VJ. The inclusion Rj C Lj is 
obvious because L/ is convex and i?/ is defined as a convex hull of a subset of L/. To 
prove the inverse inclusion, let us consider the path-connected components of D \ Rj. 
Sets Vj (j £ I) are the path-connected components of D \ Ri because they are the 
path-connected components of D \ Q, Q C Ri and Ri D Vj =0 for j £ I. There exist 
no other path-connected components of Q C Ri because all the vertices of Vi (i G I) 
belong to Ri by construction, hence, Do\Ri C UjgjVj. Due to Lemma [2TT1 every 
path-connected component of D C Ri includes at least one vertex of D. Therefore, 
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Vj (j £ I) are all the path-connected components of D \ Rj and D \ Rj = UjgjVj. 
Finally, Rj = D \ Uj&Vj = QU (U ieI Vi) = Lj. □ 

According to Lemma 12. 6\ each path-connected component Wi C D\U can be 
represented in the form Wi = Vi \U , where ^ is a path-connected component of 
D\Q. By construction, Q (ZU, hence 

Wi = (QUVi)\C/. (2.5) 

If J7 is given by a system of inequalities then representations (|2.4p and (|2.5p give us 
the possibility to represent W{ by inequalities. Indeed, the convex envelope of a finite 
set in (j2.4p may be represented by a system of linear inequalities. If the sets Q U Vi 
and U in (|2.5p are represented by inequalities then the difference between them is also 
represented by the system of inequalities. 

The description of the path-connected component of D \ U may be constructed 
by the following steps: 

1. Construct the graph of the 1-skeleton of D, this is D\] 

2. Find the vertices of D that belong to U , this is the set /To; 

3. Find the edges^of D that intersect U , this is the set U\. 

4. Delete jfrom D\ all the vertices from Uq and the edges from J7i, this is the 
graph D 1 \U; ^ 

5. Find all the connected components of D\ \ U. Let the sets of vertices of these 
connected components be Vbi, . . . , Vo q ; 

6. Select the edges e of D which intersect U but the intersection e H U does not 
include vertices of D. For every such an edge select one point x e G e HU. 
The set of these points is Q\. 

7. For every i = 1, . . . , q describe the polyhedron Ri = conv(Z7o U Q± U Vbi); 

8. There exists q path-connected components of D \ U: Wi = Ri\U. 

Every step can be performed by known algorithms including algorithms for the so- 
lution of the double description problem [49j [I4j [21] and the convex hull algorithms 

1551 . 

Let us use the simple system of three reagents, ^-1,2,3 (Fig- H-ip to illustrate the 
main steps of the construction of the path-connected components. The polyhedron D 
is here the 2D simplex (Fig. II. Ik ). The plane AffD is given by the balance equation 
C1+C2+C3 = 1. We select U = {c | G(c) < go} as an example of a convex set (Fig. 12. 2k ). 
It includes no vertices of D, hence, Uq = 0. U intersects eachjsdge of D in the middle 
point, hence, U\ includes all the edges of D. The graph Di\U consists of three 
isolated vertices. Its connected components are these isolated vertices. Qi consists of 
three points, the middles of the edges (1/2,1/2,0), (1/2,0,1/2) and (0,1/2,1/2) (in 
this example, the choice of these points is unambiguous, Fig. 12.2k ). 

The polyhedron Q is a convex hull of these three points, that is the triangle given 
in Aff(D) by the system of three inequalities ci^,3 < 1/2 (Fig. 12.2b ). The connected 
components of D \ Q are the triangles Vi given in D by the inequalities q > 1/2. In 
the whole R 3 , these sets are given by the systems of an equation and inequalities: 

Vi = {c I Ci, 2 ,3 > 0, Ci + c 2 + c 3 = 1, Ci > 1/2} . 

The polyhedron Ri is the convex hull of four points, the middles of the edges and 
the zth vertex (Fig. 12.2b ). In D, Ri is given by two linear inequalities, Cj < 1/2, j ^ i. 
In the whole R 3 , these inequalities should be supplemented by the equation and 
inequalities that describe D: 



Ri = {C I Ci, 2 ,3 > 0, Ci + C 2 + C 3 = 1, Cj < 1/2 (j ^ i)} . 
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Fig. 2.2. Construction of the path- connected components Wi of D\U for the simple example, 
(a) The balance simplex D, the set U and the path- connected components Wi; (b) The polyhedron 
Q = conv(t/o U Qi) \2. 3\) (Uo = 0, Qi consists of the middles of the edges) and the connected 
components Vi of D\ Q: Vi = {ce D\a> 1/2}; (c) The set R x = conv(U UQiU (D n Vi)); (d) 
The connected components W\ described by the inequalities (as R\\U \2.4\j )- 

The path-connected components of D \ U, Wi are described as Ri \ U Fig. 12.21 1): 
in D we get Wi = {c \ cj < 1/2 (j ^ z), G(c) > g^}. In the whole R 3 , 

Wi = {c I Ci, 2 ,3 > 0, ci + c 2 + c 3 = 1, < 1/2 (j ^ z), G(c) > £0} • 

Vi are convex sets in this simple example, therefore, it is possible to simplify 
slightly the description of the components Wi and to represent them as Vi \ U: 

Wi = {c\ ci 5 2, 3 > 0, ci + c 2 + c 3 = 1, Ci > 1/2, G(c) > g } 

(or W, = {c I q > 1/2, G(c) > g } in D). 

In the general case (more components and balance conditions), the connected 
components Vi may be non-convex, hence, description of these sets by the systems 
of linear equations and inequalities may be impossible. Nevertheless, there exists 
another version of the description of Wi where a smaller polyhedron is used instead 
of Ri _ 

Let Voi be the set of vertices of a connected component of the graph Di \ U. Let 
E ou t(Voi) be the set of the outer edges of Vo% in D\ i.e., this is the set of edges of D\ 
that connect vertices from Voi with vertices from Do \ Voi- For each e G E out (Voi) 
the corresponding edge e C D\ intersects U because Vo% is the set of vertices of a 
connected component of the graph D\\U . 

Let us select a point x e G UC\e for each e G E out (Voi) (we use the same notations 
for the edges from D\ and the corresponding edges from D\). Let us use the notation 

Qoi = {x e \e e E out (V 0i )} . 

Proposition 2.9. The path- connected component Wi of D\U allows the follow- 
ing description: 



Wi = coiiv(Q iUVoi)\U. 
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Proof. The^set Qo* C U and Voi is the set of vertices of a connected component 
of the graph Di \ Qoi by construction because Qoi cuts all the outer edges of Voi in 
D\. The rest of the proof follows the proofs of Lemma [2771 and Proposition 12.81 □ 

This proposition allows us to describe Wi by the system of inequalities. For 
this purpose, we have to use a convex hull algorithm and describe the convex hull 
conv(Qoi U Voi) by the system of linear inequalities and then add the inequality that 
describes the set \U. 

In the simple system (Fig. 12. 2p , the connected components of the graph Di\U 
are the isolated vertices. The set Qoi f° r the vertex Ai consists of two middles of its 
incident edges. In Fig. 12.2b . the set conv(Qoi U Voi) for Voi — {Ai} is the triangle V\. 

3. Thermodynamic tree. 

3.1. Problem statement. Let a real continuous function G be given in the 
convex bounded polyhedron D cW 1 . We assume that G is strictly convex in D, i.e. 
the set (the epigraph of G) 

epi(G) = {(x, g) \ x G D, g > G(x)} c D x (-oo, oo) 

is convex and for any segment [x,y] C D (x ^ y) G is not constant on [x,y]. A 
strictly convex function on a bounded convex set has a unique minimizer. Let x* be the 
minimizer of G in D and let g* = G(x*) be the corresponding minimal value. The level 
set S g = {x e D \ G(x) = g} is closed and the sublevel set U g — {x G D \ G(x) < g} 
is open in D. The sets S g and D\U g are compact and S g C D\U g . 

Let x : y G D. According to Corollary 13.41 proven in the next subsection, an 
admissible path from x to y in D exists if and only if ir(y) belongs to the ordered 
segment [tt(x*) : tt(x)]. Therefore, to describe constructively the relation x ^ y in D 
we have to solve the following problems: 

1. How to construct the thermodynamic tree T? 

2. How to find an image ir(x) of a state x G D on the thermodynamic tree T? 

3. How to describe by inequalities a preimage of an ordered segment of the 
thermodynamic tree, 7r -1 ([iu, z\) C D (w, z G T, z ^ w)l 

3.2. Coordinates on the thermodynamic tree . We get the following lemma 
directly from Definition ll.il Let x,y G D. 

Lemma 3.1. x ~ y if and only if G(x) = G(y) and x and y belong to the same 
path- connected component of S g with g = G(x). 

The path-connected components of D \ U g can be enumerated by the connected 
components of the graph D\\U g . The following lemma allows us to apply this result 
to the path-connected components of S g . 

Lemma 3.2. Let g > g* , W\, . . . , W q be the path- connected components of D\U g 
and let cr l5 . . . , a p be the path- connected components of S g . Then q = p and <r^ may 
be enumerated in such a way that ai is the border ofWi in D. 

Proof. G is continuous in D, hence, if G{x) > g then there exists a vicinity of x 
in D where G(x) > g. Therefore G{y) = g for every boundary point y of D \ U g in D 
and S g is the boundary of D \ U g in D. 

Let us define a projection g : D \ U g —> S g by the conditions: O g {x) G [x,x*] 
and G(6 g (x)) — g. By definition, the inequality G{x) > g holds in D \U g . The 
function f x (X) = G((l — A)x* + Ax) is strictly increasing, continuous and convex 
function of A G [0,1], f x (0) = g* < g, f x (l) = G(x) > g. The function f x (X) 
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depends continuously on x G D \ U g in the uniform metrics. Therefore, the solution 
\ x to the equation f x (X) = g on [0,1] exists (the intermediate value theorem), is 
unique, and continuously depends on x G D \U g . The projection 6 g is defined as 
g (x) = (1 - X x )x* +\ x x. 

The fixed points of the projection 6 g are elements of S g . The image of each 
path-connected component Wi is a path-connected set. The preimage of every path- 
connected component G{ is also a path-connected set. Indeed, let g (x) G <ji and 
g (y) G <Ji. There exists a continuous path from x to y in D\ U g . It may be composed 
from three paths: (i) from x to g (x) along the line segment [x, 6 g {x)\ C [x, x*] 
then a continuous path in (Ji between g (x) and g (y) (it exists because &i is a path- 
connected component of S g and it belongs to D\U g because S g C D\U g ) and, finally, 
from g (y) to y along the line segment [0 g (y),y] C [x*,y]. Therefore, the image of 
a path-connected component Wi is a path-connected components of S g that may be 
enumerated by the same index z, G{. This &i is the border of Wi in D. □ 

The equivalence class of x G D is defined as [x] = {y G D\y ~ x}. Let W(x) be 
a path-connected component of D \ U g (g = G{x)) for which 6 g (W(x)) = [x]. Due to 
Lemma [321 such a component exists and 

W(x) = {yeD\y>zx}. (3.1) 

Let us define a one-dimensional continuum y that consists of the pairs (g,M), 
where g* < g < ^ max and M is a set of vertices of a connected component of D\\U g . 
For each (g, M) the fundamental system of neighborhoods consists of the sets V p 
(P>0): 

V p = {(</, MO | (^, MO ey,\g- g'\ < p, M' C M} . (3.2) 
Let us define the partial order on 3^: 

(flS M ) (9', M ') '^9>9 f and M C M' . 
Let us introduce the mapping uj : D — >• 3^: 

w(a?) = (G(x), W(x) HL>o). 

Theorem 3.3. T/iere exzs^ a homeomorphism between y and T that preserves 
the partial order and makes the following diagram commutative: 

UJ 

y 

Proof. According to Lemmas 13. 2\ 13.11 and Proposition \2A\ uo maps the equivalent 
points x to the same pair (g,M) and the non-equivalent points to different pairs 
(<7, M). For any y G D, x £3 y if and only if cj(x) £3 a; (2/). 

The fundamental system of neighborhoods in Y may be defined using this partial 
order. Let us say that (g, M) is compatible to (g' , M') if (g' , M') £3 (g, M) or (g, M) >^ 
(g',M')}. Then for p > 




F p = {(^, MO ey\\j- j'\ < p and (#', MO is compatible to (g, M)} . 
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For sufficiently small p this definition coincides with (|3.2p . 

So, by the definition of T as a quotient space D/ ~, y has the same partial 
order and topology as T. The isomorphism between y and T establishes one-to- 
one correspondence between the tt- image of the equivalence class [x], tt([x]), and the 
cj-image of the same class, u; ([#]). □ 

y can be considered as a coordinate system on T. Each point is presented as a 
pair (g, M) where g* < g < g ma x and M is a set of vertices of a connected component 
of Di\U g . The map uj is the coordinate representation of the canonical projection 
7r : D — > T. Now, let us use this coordinate system and the proof of Theorem 13.31 to 
obtain the following corollary. 

Corollary 3.4. An admissible path from x to y in D exists if and only if 

n(y) e [7r(x*),7r(x)] . 

Proof Let there exist an admissible path from x to y in D, cp : [0, 1] — >• D. 
Then tt(x) ^ n(y) in T. Let 7r(x) = (G(x),M) in coordinates y. For any u G M, 
7r(y) G [7r(x*),7r(v)] and 7r(x) G [7r(x*), 7r(v)j. 

Assume now that n(y) G [7r(x*),7r(x)] and7r(x) = (G(x),M). Then the admissible 
path from x to y in D can be constructed as follows. Let v G M be a vertex of D. 
G(i>) > G(x) for each v G M. The straight line segment [x*,v] includes a point x\ 
with G(x\) = G(x) and i/i with G(yi) = G(y). Coordinates of tt(xi) and tt(x) in y 
coincide as well as coordinates of 7t(t/i) and ir(y). Therefore, x ~ x\ and y ~ y 1 . The 
admissible path from x to y in D can be constructed as a sequence of three paths: 
first, a continuous path from x to x\ inside the path-connected component of Sq{x) 
(Lemma 13. ip . then from x\ to y\ along a straight line and after that a continuous 
path from y\ to y inside the path-connected component of Sg( v )- D 

To describe the space T in coordinate representation y, it is necessary to find the 
connected components of the graph D\ \ U g for each g. First of all, this function, 

g ^ the set of connected components of D\ \ U g , 

is piecewise constant. Secondly, we do not need to solve at each point the computa- 
tionally heavy problem of the construction of the connected components of the graph 
D\ \ U g "from scratch" . The problem of the parametric analysis of these components 
as functions of g appears to be much cheaper. Let us present a solution of this prob- 
lem. At the same time, this is a method for the construction of the thermodynamic 
tree in coordinates (g,M). 

The coordinate system y allows us to describe the tree structure of the continuum 
T. This structure includes a root, (^*,.Do), edges, branching points and leaves. 

Let M be a connected component of Di\U g for some g, g* < g < g max - If M ^ Do 
then the set of all points (g,M) G T has for a given M the form (a M ,%] x M . 
Qlm < ^m- We call this set an edge of T. 

If M includes all the vertices of D (M = Do) then the set of all points (g, M) G T 
has the form [g*,aD ] x ^o- This may be either an edge (if cld > g*) or just a root, 
{(g*,Do)}, (this is possible in ID systems). 

Let us define the numbers a M = inf{g | (g, M) G T}. Let us introduce the set of 
outer edges of Min Di, E out (M). This is the set of edges of Di that connect vertices 
from M with vertices from Dq \ M . We keep the same notation, E out (M), for the set 
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of the corresponding edges of D. 

a M — max min{G(x) | x G e} . (3.3) 

eeE out (M) 

This number, a M , is the "cutting value" of G for M. It cuts M from the other vertices 
of Di in the following sense: if we delete from D\ all the edges e with the label values 
< a M then M will remain attached to some vertices from D \M. If we delete the 
edges with the label values < a M then M becomes disconnected from Do\M. There 
is the only connected component of D\ \ Ua M that includes M, M' ^ M. The pair 
(a M , M') G T is a branching point of T. The edge (a M ,%] x M connects two vertices, 
the upper vertex (aM,M) and the lower vertex, (a M ,M / ). 

If M consists of one vertex, M = {v}, then the point (G(v), {v}) is a leaf of T. 

3.3. Construction of the thermodynamic tree. To construct the tree of G in 
D we need the graph D\ of the 1-skeleton of the polyhedron D. Elements of D\ should 
be labeled by the values of G. Each vertex v is labeled by the value j v = G(v) and 
each edge e = [v,w] is labeled by the minimal value of G on the segment [v,w] C D, 
g e = min[ v w ] G(x). We need also the minimal value g* = min^{G(x)} because the 
root of the tree is (g*,Do). 

The strictly convex function G achieves its local maxima in D only in vertices. 
The vertex v is a (local) maximizer of g if g e < j v for each edge e that includes v. 
The leaves of the thermodynamic tree are pairs (7^, {v}) for the vertices that are the 
local maximizers of G. 

As a preliminary^ step of the construction, we arrange and enumerate the labels 
of the elements of Di, the vertices and edges, in descending order. Let there exist I 
different label values: g max = a\ > a2 > . . . > a/. Each a& is a value j v = G(v) at 
a vertex v G Z^o or the minimum of G on an edge ecDi (or both). Let A{ be the 
set of vertices v G Do with j v = ai and let E{ be the set of edges of D\ with g e = ai 

c 1 /)• _ 

Let us construct the connected components of the graph D\ \ U g starting from 
&i = ^max- The function G is strictly convex, hence, a\ = j v for a set of vertices 
A\ C Do but it is impossible that a\ = ^e_for an edge e, hence, E\ = 0. 

The set of connected components of D\ \ U g is the same for all g G (a^+i, a*]. For 
an interval (0,2, cti] the connected components of D\ \ U g are the one-element sets {v} 
for v e Ai. 

For ^ G [<7*,a/] the graph \ U g includes all the vertices and edges of D\ and, 
hence, it is connected for this segment. Let us take, formally, a^ + i = g*. 

Let Ci = {M|, . . . , M^,.} be the set of the connected components of D\ \ U g for 
g G (a^, di-i] (i = 1, . . . , I). Each connected component is represented by the set of 
its vertices Mj. Let us describe the recursive procedure for construction of Ci'. 

1. Let us take formally Co = 0. 

2. Assume that C{-\ is given and i < I. Let us find the set Ci of connected 
components of D\ \ U g for g = (and, therefore, for g G (a^+i, a^]). 

• Add the one-element sets {v} for all ^ G A« to the set 

^^{Mr 1 ,...^-^}. 

Denote this auxiliary set of sets as L^o = {Mi,...,M g }, where g = 
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• Enumerate the edges from E{ in an arbitrary order: ei, . . . , ei^j . For 
each = 0, . . . , \Ei\, create recursively an auxiliary set of sets L^k by the 
union of some of elements of Lk-i'. Let L^k-i be given and e& connects 
the vertices u and If v and belong to the same element of L^-i 
then = L^k-i- If ^ and belong to the different elements of L^k-i, 
M and M', then L iy ^ is produced from L^k-i by the union of M and 
M'\ 

L z ,k = (L <jfc _i \ {M} \ {M'» U {M U M'} 

(we delete two elements, M and M', from L^k-i and add a new element 
MUM'). 

The set A of connected components of D\ \ U g for g = en is Ci = Li^ E .\. 
Generically, all the labels of the graph Di vertices and edges are different and the sets 
Ei and Ai include not more than one element. Moreover, for each i either E{ or Ai 
is generically empty and the description of the recursive procedure may be simplified 
for the generic case: 

1. Let us take formally Co = 0. 

2. Assume that Ci-\ is given and i < I. 

• If ai is a label of a vertex v, ai = j v , then add the one-element set {v} 
to the set Ci-\\ Ci = Ci-\ U {{^}}. 

• Let ai be a label of an edge e = [v,?/]. If v and v' belong to the same 
element of A_i then Ci = A-i- If v an d v' belong to the different 
elements of M and M r , then Ci is produced from Ci-\ by the 
union of M and M' (delete elements M and M' and add an element 
MUM'). 

The described procedure gives us the sets of connected components of D\ \ U g 
for all g and, therefore, we get the tree T. The descent from the higher values of G 
allows us to avoid the solution of the computationally more expensive problem of the 
calculation of the connected components of a graph at any level of G. 

3.4. The problem of attainable sets. In this section, we demonstrate how to 
solve the problem of attainable sets. For given x G D (an initial state) we describe 
the attainable set 

Att(x) = {yeD\x>zy} 

by a system of inequalities. Let the tree T of G in D be given and let all the pairs 
(g, M) G T be described. We also use the notation Att(z) for sets attainable in T 
from z G T. 

First of all, let us describe the preimage of a point (g, M) G T in D. It can be 
described by the equation G(x) = g and a set of linear inequalities. For each edge e we 
select a minimizer of G on e, x e = argmin{G(x) | x G e} (we use the same notations 
for the elements of the graph D\ and of the continuum D\). Let 

Qm = {x e I e G E out (M)} . 

In particular, a M = max{G(x) | x G Qm}- 

The following Proposition is a direct consequence of Proposition 12.91 
Proposition 3.5. The preimage of (g,M) in D is a set 



7T" 1 (g, M) = {x G conv(Q M U M) \ G(x) = g} . 



(3.4) 
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bs [S] b s [S] b s [S] 



Fig. 4.1. The balance polygon D on the plane with coordinates [S] and [ES] for the four- 
component enzyme-substrate system S, E, ES P with two balance conditions, bs = [S] + [ES] + [P] = 
const and bs = [E] + [ES] = const. 



The sets M and Qm in (|3.4p do not depend on the specific value of g. It is 
sufficient that the point (g, M) G T exists. 

Let us considerable second projection of T, i.e., the set of all connected compo- 
nents of the graph D\ \ U g for all g. For a connected component M, the lower chain 
of connected components is a sequence M = M\ ^ M2 ^ . . . ^ M/c . ( "Lower" here 
means the descent in the natural order in T, For a given initial element M = Mi 
the maximal lower chain of M is the lower chain of M that cannot be extended 
by adding new elements. By construction of connected components, the maximal 
lower chain of M is unique for each initial element M. In the maximal lower chain 

For each set of values H C (cl m ,cim} the preimage of the set H x M C T is given 
by (pT4j) as 

tt-^H xM) = {xe conv(Q M U M) \ G(x) e H} . (3.5) 

We describe the set Att(x) for x G D by the following procedures: (i) find the 
projection tt(x) of x onto T, (h) find the attainable set in T from tt(x), AU(tt(x)), 
and (iii) find the preimage of this set in D: 

Att(x) = 7T' 1 (AU(tt(x))) . (3.6) 

The attainable set Att(g, M) in T from (g, M) G T is constructed as a union of 
edges and its parts. Let M = Mi ^ M2 § ■ • • § Affe = £>o be the maximal lower 
chain of M. Then 

AU(g,M)=(a 1 ,g] xM 1 U{a 2 ,a 1 ] x M 2 U . . . 

U (afc_i,afc_ 2 ] x Afjfe-i U [a fc ,a fc _i] x M fe , 

where a i = a M . . 

To find the preimage of Att(g, M) in D we have to apply formula (|3.5p to each 
term of (j3.7p . In Sec. 11.31 we demonstrated how to find tt(x). Therefore, each step of 
the solution of the problem of attainable set (j3.6ft is presented. 

4. Chemical thermodynamics: examples. 

4.1. Skeletons of the balance polyhedra. In chemical thermodynamics and 
kinetics, the variable Ni is the amount of the ith component in the system. The 
balance polyhedron D is described by the positivity conditions Ni > and the balance 
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Fig. 4.2. The graph D\(b) of the one-skeleton of the balance polyhedron for the six-component 
system, H2, O2, H ; O, H^O, OH ; as a piece-wise constant function of b = (bn,bo). For each vertex 
the components are indicated which have non-zero concentrations at this vertex. 



conditions b{(N) = const (i = l,...,m). Under the isochoric (the constant 

volume) conditions, the concentrations q also satisfy the balance conditions and we 
can construct the balance polyhedron for concentrations. Sometimes, the balance 
polyhedron is called the reaction simplex with some abuse of language because it is 
not obligatory a simplex when the number m of the independent balance conditions 
is greater than one. 

The graph D\ depends on the values of the balance functionals b{ = bi(N) = 
Ej=i a i^j- For the positive vectors N, the vectors b with coordinates b{ = b{(N) 

form a convex polyhedral cone in R m . Let us denote this cone by A. D\(b) is a piece- 
wise constant function on A. Sets with various constant values of this function are 
cones. They^ form a partition of A. Analysis of this partition and the corresponding 
values of D\ can be done by the tools of linear programming [26] . Let us represent 
several examples. 

In the first example, the reaction system consists of four components: the sub- 
strate S, the enzyme E, the enzyme-substrate complex ES and the product P. we 
consider the system under constant volume. We denote the concentrations by [5], 
[E], [ES] and [P]. There are two balance conditions: b s = [S] + [ES] + [P] = const 
and b E = [E] + [ES] = const. 

For bs > b E the polyhedron (here the polygon) D is a trapezium (Fig. 14.1k ). 
Each vertex corresponds to two components that have non-zero concentrations in this 
vertex. For bs > bE there are four such pairs, (ES,P), (ES,S), (E,P) and (E,S). 
For two pairs there are no vertices: for (S,P) the value bE is zero and for (ES,E) it 
should be bs < bE- When bs = &e, two vertices, (ES, P) and (ES, 5), transform into 
one vertex with one non-zero component, ES, an the polygon D becomes a triangle 
(Fig. 14.1b ). When bs < bE then D is also a triangle and a vertex ES transforms in 
this case into (ES,E) (Fig. Oh). 
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Fig. 4.3. Transformations of the graph D±(b) with changes of the relation between 6h and bo: 
(a) transition from the regular case 6h > 26o to the regular case 2bo > 6h > fro through the singular 
case 6h — 26o ; transition from the regular case 2bo > bn > bo to the regular case 6h < fro 
through the singular case 6h = bo- 



For the second example, we select a system with six components and two balance 
conditions: H 2 , 2 , H, O, H 2 0, OH; 

b u = 27V H2 + N u + 27V H2 o + ^oh , 
b = 2No 2 + No + 7Vh 2 o + A^oh . 

The^cone A is a positive quadrant on the plane with the coordinates 6h, &o- The 
graph Di(b) is constant in the following cones in A (6h,^o > 0): (a) 6h > 2bo, (b) 
6 H = 26 , (c) 2b >b u >b , (d) 6 H = b and (e) 6 H < b (Fig. S2J). 

The cases (a) 6 H > 26 , (c) 2b > bn > b , and (e) 6 H < (Fig. [472]) are 
regular: there are two independent balance conditions and for each vertex there are 
exactly two components with non-zero concentration. In case (a) (&h > 2&o)> if 
6h — >• 2bo then two regular vertices, H 2 , H 2 and H, H 2 0, join in one vertex (case 
(b)) with only one non-zero concentration, H 2 (Fig. 14.3k ). This vertex explodes in 
three vertices O, H 2 0; 2 , H 2 and H 2 0, OH, when 6h becomes smaller than 2bo 
(case (c), 2bo > > bo) (Fig. 14.3k ). Analogously, in the transition from the regular 
case (c) to the regular case (e) through the singular case (d) (6r = 60) three vertices 
join in one, OH that explodes in two (Fig. 14.3b ). 

For the modeling of hydrogen combustion, the eight-component_model is used 
usually: H 2 , 2 , H, O, H 2 0, OH, H 2 2 , H0 2 . In Fig.0~4]tlie graph D 1 is presented 
for one particular relations between 6h and 26o, &h = 2&o- This is the so-called 
"stoichiometric mixture" where proportion between bu and 2bo is the same as in the 
"product", H 2 0. 
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Fig. 4.4. The graph D± for the eight- component system, H 2 , O2, H, O, H 2 ; OH 7 H 2 02 7 HO2 
for the stoichiometric mixture, 6h = 26o- T/ie vertices that correspond also to the six-component 
mixture are distinguished by bold font. 



4.2. Examples of the thermodynamic tree. In this section, we present two 
example of the thermodynamic tree. First, let us consider the trapezium (Fig. 14.1k ). 
Let us select the order of numbers j v and g e according to Fig. 14.51 The vertices and 
edges are enumerated in order of j v and g e (starting from the greatest values). The 
tree is presented in Fig. 14.51 On the right, the graphs D\ \ U g are depicted for all 
intervals (a^_i,a^]. For (72,71] it is just a vertex v\. For (#1,72] it consists of two 
disjoint vertices, v\ and vi> For (73,01] these two vertices are connected by an edge. 
On the interval (#2? 73] the graph D± \ U g is an edge (^1,^2) and an isolated vertex 
vs. On (74,^2] all three vertices v±, V2 and ^3 are connected by edges. For (03,74] 
the isolated vertex V4 is added to the graph D\ \ U g . For g < gs the graph D\ \ U g 
includes all the vertices and is connected. 

For the second example (Fig. 14.6)) we selected the six-component system (Fig. 14.2)) 
with the stoichiometric hydrogen-oxygen ratio, 6h = 26o- The selected order of 
numbers 7$, gj is presented in Fig. 14.61 

5. Conclusion. We studied dynamical systems that obey a continuous strictly 
convex Lyapunov function G in a positively invariant convex polyhedron D. Convexity 
allows us to transform n-dimensional problems about attainability and attainable sets 
into an analysis of ID continua and discrete objects. 

We construct the tree (the Adelson-Velskii - Kronrod - Reeb tree [1, 42, 56]) of 
the function G in D and call this ID continuum the thermodynamic tree. 

The thermodynamic tree is a tool to solve the "attainability problem" : is there a 
continuous path between two states, x and y along which the conservation laws hold, 
the concentrations remain non-negative and the relevant thermodynamic potential G 
(Gibbs energy, for example) monotonically decreases? This question arises often in 
non-equilibrium thermodynamics and kinetics. The analysis of the admissible paths 
can be considered as a dynamical analogue of the study of the steady states feasibility 
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E, P e A E,S 

7l>72>gl > 73 > g2 > 74 > g3 > g4 




ES,P 



E,P 



Fig. 4.5. The thermodynamic tree for the four- component enzyme-substrate system S, E, ES 
P (Fig. |^.ip with excess of substrate: bs > bs (case (a)). The vertices and edges are enumerated 
in order of j v and g e (starting from the greatest values). The order of these numbers is represented 
in Fig. On the right, the graphs D\ \ U g are depicted. The solid bold line on the tree is the thermo- 
dynamically admissible path from the initial state E,S (enzyme plus substrate) to the equilibrium. 
There are leaves at all levels g = ji- There are branching points at g = #1,2,3 o/nd no vertices at 
9 = 94. 



in chemical and biochemical kinetics. In this recent study, the energy balance method, 
the stoichiometric network theory, the entropy production analysis and the advanced 
algorithms of convex geometry of polyhedral cones are used [5J [54] . 

The obvious inequality, G(x) > G(y) is necessary but not sufficient condition for 
existence of an admissible path from x to y. In ID systems, the space of states is an 
interval and the thermodynamic tree has two leaves (the ends of the interval) and one 
root (the equilibrium). In such a system, a spontaneous transition from a state x to 
a state y is allowed by thermodynamics if G(x) > G(y) and x and y are on the same 
side of the equilibrium, i.e. they belong to the same branch of the thermodynamic 
tree. This is just a well known rule: "it is impossible to overstep the equilibrium in 
one-dimensional systems" . 

The construction of the thermodynamic tree gives us the multidimensional ana- 
logue of this rule. Let 7r : D — >> T be the natural projection of the balance polyhedron 
D on the thermodynamic tree T. A spontaneous transition from a state x to a state 
y is allowed by thermodynamics if and only if 7r(y) G [^(x), 7r(iV*)], where N* is the 
equilibrium and [ir(x),ir(N*)] is the ordered segment. 

In this paper, we developed methods for solving the following problems: 

1. How to construct the thermodynamic tree T? 

2. How to solve the attainability problem? 

3. How to describe the set of all states attainable from a given initial state xl 
For this purpose, we analyzed the cutting of a convex polyhedron by a convex 

set and developed the algorithm for construction of the tree of level set components 
of a convex function in a convex polyhedron. In this algorithm, the restriction of G 
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Fig. 4.6. T/ie thermodynamic tree for the six-component H2-O2 system, H2, O2, H 7 O, H20 ; 
OH with the stoichiometric hydrogen-oxygen ratio, bn = 26q (Fig. \4~!fy )- The order of numbers ^i, 
gj is presented in Fig. On the right, the graph D\ \ U g is represented for g = gio. For g < gio, 
the graph D\ \ U g includes all the vertices and is connected. The solid bold line on the tree is the 
thermodynamically admissible path from the initial state H2, O2 to the equilibrium. There are leaves 
at all levels g = ji- There are branching points at g = gi— 4,7,10 and no vertices at g = #5,6,8,9- 



onto the 1-skeleton of D is used. This finite family of convex functions of one variable 
includes all necessary information for analysis of the tree of the level set component 
of the convex function G of many variables. 

In high dimensions, some steps of our analysis become computationally expensive. 
The most expensive operations are the convex hull (description of the convex hull of 
a finite set by linear inequalities) and the double description operations (description 
of the faces of a polyhedron given by a set of linear inequalities). Therefore, in high 
dimensions some of the problem may be modified, for example, instead of the explicit 
description of the convex hull it is possible to use the algorithm for solution of a prob- 
lem: does a point belong to this convex hull [55]. The computational aspects of the 
discussed problems in higher dimensions deserve more attention and the proper mod- 
ifications of the problems should be elaborated. For example, two following problems 
need to be solved efficiently: 

• To find the maximal and the minimal value of any linear function / in a class 
of thermodynamic equivalence; 

• To evaluate the maximum and the minimum of dG/dt in any class of ther- 
modynamic equivalence: —a < dG/dt < —a < 0. 

For any w G T, the solution of the first problem allows us to find an interval of 
values of any linear function of state in the corresponding class of thermodynamic 
equivalence. We can use the results of Sec. 12.21 to reformulate this problem as the 
convex programming problem. 

The second problem gives us the possibility to consider dynamics of relaxation 
on T. On each interval on T we can write 



-a(g) <dg/dt < -a(g) <0, 



(5.1) 
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where the functions cr(g)^a(g) > depend on the interval on T. 

This differential inequality ([5.1)) will be a tool for the study of the dynamics of re- 
laxation and may be considered as a reduced kinetic model that substitutes dynamics 
on the d-dimensional balance polyhedron D by dynamics on the one-dimensional den- 
drite. The problem of the construction of the reduced model ([5.1)) is closely related to 
the following problem [66 : "Can one establish a lower bound on the entropy produc- 
tion, in terms of how much the distribution function departs from thermodynamical 
equilibrium?" In 1982, C. Cercignani [12 proposed a simple linear estimate for a(g) 
for the Boltzmann equation (Cercignani's conjecture). After that, these estimates 
were studied and improved by many authors [15j[ll ; 65 , 66 and now the state of art 
achieved for the Boltzmann equation gives us some hints how to create the relaxation 
model (|5.1|) on the thermodynamic tree for the general kinetic systems. This may be 
the next step in the study of the thermodynamic trees. 
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